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Abstract 

We discuss a model for interferometric GW antennas illuminated by a laser beam and 
a vacuum squeezed field. The sensitivity of the antenna will depend on the properties of 
the radiation entering the two ports and on the optical characteristics of the interferometer 
components, e.g. mirrors, beam-splitter, lenses. 


1 Introduction 

An important ingredient for improving the sensitivity of Michelson interferometric gravitational 
wave detectors (GWD) is using appropriate states for the light beams illuminating its two input 
ports. In interferometric measurements the quantum noise is due to the fluctuations of the number 
of photons and to the random motion of the mirrors induced by the radiation pressure. The GW 
signal is extracted from the spectral density of the output. 

Purpose of this paper is to discuss the dependence of the sensitivity of an interferometric 
GW antenna on the photon-noise and radiation pressure noises. In particular we will consider 
an interferometer driven by a fluctuating laser beam and a squeezed-vacuum field generated by a 
degenerate 0P0 driven by the second harmonic of the laser beam. Particular attention will be 
paid to the influence of phase and amplitude fluctuations of the laser beam. 


2 Michelson interferometer 

We consider a Michelson interferometer with two mirrors Mi and M 2 suspended at the ends of 
two arms. The vertices of Mi and M 2 are located on the axes y and x passing through the origin 
0, while the beam splitter is centered on 0 (see Fig. 1). 

In order to account for aberration effects, we will model the interferometer as a multimode 
device: we consider two groups of beams entering through the ports Pi, P 2 (Fig. 1), described by 

1 also with Istituto Nazionale di Ottica (INO) 
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Figure 1: Michelson’s interferometer. B.S.= beam-splitter; Wm aberration regions due to the 
mirrors; Wb.s. aberration regions due to the beam-splitter; F.D.= frequency doubler; O.P.O.= 
Optical Parametric Oscillator; di,,^, = in-fields at the port 1,2; bi t ,b 2 , = fields at the mirror 1,2; 
bi,,c 2 , = out-fields at the port 1,2. 


the operators (dj, , aj, ), (d 22 , dj 2 ) with j = 0, 1, ..., N-l, acting on a Hilbert space K a = K ai ®K a2 , 
with K ai = K^ <g> . . . ® K[^~ l \ * — 1,2. More specifically, the modes relative to port P\ consist 
of Gauss-Hermite beams travelling along the x-axis with waist in O, 


u lm (y,z-,x) oc e-^H, ff„ 



The pair of indices Im will be denoted by ij. Analogously, for P 2 we consider a similar family of 
Gauss-Hermite beams propagating along the y-axis with waist in O (see Fig. 2). 

Passing through the beam-splitter, the input beams transform in two fields at M \ , M 2 described 

by {bji »&],)» (hi > ) with J - 0, 2, . . . , N - 1 , acting on a Hilbert space Kb = Kb x ® K ^ , such that 

Kb, = K { £ ] ®. . .®K[^~^ and two outgoing beams described at Pi, P 2 by (dj, ,cj,), (c J2 , cj 2 ) with j = 
0, 1, . . . , N — 1, acting on a Hilbert space K c = K Cl ®K C2 , where, again, K c , = K$ ® . . . ® K[^~ l \ 

For the sake of notational convenience, we introduce the bold symbols a, = d u ® a, 2 = 

for indicating the pair of i-th modes relative to Pi and P 2 respectively. Analogously, we introduce 
/ d 0 \ 

the vector A = • . With the same meaning, we will introduce the vectors B and C for 

^ &N - 1 / 

the operators bj and Cj, relative respectively to the mirrors and the output ports. 
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Figure 2: Schematic of Gauss-Hermite beams. 


Assuming the fraction of energy lost during the passage through the interferometer be indipen- 
dent of the mode considered, B and C can be redefined as vectors proportional to the actual ones 
and carrying the same energy of A. In view of the energy conservation, the linearity and time 
invariance of the antenna, the outgoing vector C can be related to the ingoing one A by the 
unitary matrix U, 

C = U-A, (1) 


with 


U = 


( Uo,o tfo.1 

0i.o 0i.i 

l 0jv - 1,0 0JV-U 


00, JV— 1 \ 

01. JV-1 

0Af— l.JV— 1 } 


(2) 


where each Uij is a 2 x 2 matrix. Moreover, to preserve the bosonic commutation relations, U 
must be a symplectic matrix, that is U £ Sp(2N, R) and Uij £ Sp( 2, R). 

From now on we will consider an interferometer illuminated by two T EMoo gaussian modes 
on ports P\ and P 2 respectively. This amounts to considering an input state vector of the form 


|*>=|0o>|Oi>...|Ojv-1>, (3) 

where |0, >= |0,, > |(h 2 > indicates that the modes ij, i 2 are unexcited ground state. As a result 
of the propagation through the imperfect interferometer, the states of all these modes will be 
mixed up to some extent. So that, a mode initially in the ground state will be partially excited 
at the output ports. 

In view of (3), it is worth splitting R a in the product R a = R ao ® Rap, whith R ao relative to 
the fundamental modes entering the two ports, and R a& relative to the remaining 2 N — 2 modes. 
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In the same manner we 


will write A — f ?° ] , where do = ( «°' ] , dp — 

\ a 0 J \ a <h J 


( fl! \ 


/ 


, and so for 


B and C. 

Analogously to (1), the fields B at the mirrors will depend linearly on A 


B = V ■ A, 


whith V a unitary linear transformation V 


( Vi Vop \ 

V Vio v fi J 


Physically, V describes reflection 


and transmission at the beam-splitter, followed by propagation through the interferometer arms. 
Then, it can be expressed as the product 


V = $ • $(bs) • K, 

where K describes the aberration-free beam-splitter, $ ( bs) is the aberration matrix relative to the 
beam-splitter itself and $ is the interferometer arm delay matrix. 

Introducing the N x N matrices 1 X} = &i x j, Si 3 j 3 , li \j = » an< ^ the 2x2 matrices fcj, k 2 , 

we can write K, whose elements are 2x2 matrices, as 

K = e'^icos 7 <7 0 1 + sin 7 <r x I). (4) 


with <r 0 ,<7 1 Pauli matrices. The aberrations of the beam splitter are modeled by including at 
the two output faces two transparencies characterized by the aberration eykonals W( BS ){x, z) and 
W(BS){y, z ) f° r the faces perpendicular to y- and x-axes respectively. $(bs), representing the 
aberration eykonal phase factor, is symmetric with respect to the exchange of the pair of indices 
»i,*2 with ji, j 2 . 




f $(BS)0,0 

$(BS)0,1 

$(BS)0,N-1 ^ 


$(BS) = 

$(BS) 1,0 

$(BS)1,1 

$(BS)1,N-1 


: 

j 

| • 



\ $(BS)AT-1,0 

$(BS)N-1,1 

••• $(BS)JV-l,iV-l / 

where 


${BS)ij = ( 

e * w (flS)ilJi 

0 



0 



The symmetry of the $ ( bs) matrix is a consequence of the identities VP( B S). i,.j, 

and VP(BS)i 2 ,j 2 = W(BS)j 3 ,i 3 ■ 

Finally, the time-delay matrix is diagonal 


(5) 


( 6 ) 


W(BS)]uh 


$ = 



0 


! 1 


\ 0 0 ... 1 / 


(7) 
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, 4>i t 2 representing the phase delay of the *i,2-th Gauss-Hermite mode 

hitting A/i, 2. In particular, 

2 = ^•^' 1,2 ^$>1,2 "h "~f & 4 > (aua)i '2 T T ^^(pres) 1,2 A" ^^(rp)!^ (^) 

where S4>i 13 (> 0 ) stands for the delay of the i-th Gauss-Hermite mode with respect the phase 
delay fcA 1|2 of a plane wave. Scpow 1,2 represents the gravitational wave (S(f)GW } = 4 >gw 2 )• The 

other terms stand for : (i) <$0(, usp )=noise transmitted to the mirrors through the suspensions, (ii) 
5<t>{mir)= noise caused by the vibration modes of the mirrors, (iii) S<t>( prea )= pressure fluctuations 
in the partially evacuated pipes of the interferometer arms, and (iv) 5(f)( rp )= radiation pressure 
noise. 

As a result of the reflection on M\ and M 2 , the different modes propagate toward the exit 
ports, by retracing the same paths followed before. Then, 

U = —K ■ $(bs) ■ $ • $(Af) • $ • $(bs)K ( 9 ) 

where $(Af) is the mirror aberration matrix. 

3 Interferometer output 

The interferometer output is proportional to the expectation value of the difference I between the 
photocurrents detected at the ports 1 and 2 respectively 

/ = £( C U> - cU) = ct-*3-C = Al-S-A (10) 

1 

where S is the unitary self-adjoint operator S = ■ 03 • U. 

Introducing the quantity A'3 = K * • <73 • K = cos(27)<T3 1 — sin(2~f)<72 1 (see Eq. 4 ) and 
assuming an input state of the form ( 3 ), it yelds 

S = A'* • ■ 4>(bs) • A 3 • $(£s) • • $(M) ■ $ • $<bs) • K = Sjo) + S( a b) (H) 

where 5(0) = —sin <f>cri — cos <£<73, with 4> = 2(<po 2 — <Aj,), is the matrix in absence of aberrations 
and S~f = 0, while S( a b) = c*i<7i + a 2 er 2 + a 3 < 7 3 describes the effects of the aberrations and the 
deviation from the condition of exact equipartition of the incident intensity between the two B.S. 
outputs. 

Now introducing the quantities A\ = • <T\ • Oq , A 2 = a| • <7 2 • «o 1 A3 = « 0 • <r 3 ■ ao 

and considering an interferometer operating on a dark fringe, we can express the photocurrent 1 
as / = / (< i) + Av, that is as the sum of a deterministic part 1(d) =< A 3 > [S<j>GW A £*3] depending 
on the GW signal and the aberrations, and a noise depending part 

Av ~ A[( — 1 A «i) A A 2 ot 2 ~i~ < A 3 > A (A3— < A3 >)o: 3 

= /jV (pn) + ^W) + fym.r) + ^(pr«) + ^'(rp) + ^(«6) ( 12 ) 


with*, = ('*" ° j 
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In particular, as regards to the radiation pressure noise In { tp) , the mirrors Mi and M 2 can be 
considered as multiple damped pendula driven by known time dependent pressure forces, 

yi(t),x 2 (t) = f ^ (4, 2 (OMO + r M (t - t')dt' (13) 

with the impulse response of the mirrors. Accordingly 8<t>( rp )(t ) = k(y x (t) — x 2 {t)) = 

r M * (6t • cr 3 • 6), having indicated with Tm* the convolution integral (13). So That lN {rp) = 
(— l) fc+1 (rAf * A 2 ) < A 3 > . 


4 Fourier analysis of the interferometer output 


In most GW antennas the signal is extracted from the frequency spectrum of the photocurrent 
I = 1(d) + In- Therefore, the sensitivity of the interferometer depends on the autocorrelation of I, 

< : I( T )iI(0) :>=<: 1(d) y 1(d) '■> + <: In, In •> (14) 


having considered 1(d) and In as indipendent. 

The limiting sensitivity of the antenna will be obtained by equating the Fourier component 
Sgw(u) of <: 1(d), 1(d) ■> at the frequency of the gravitational wave to the noise component, 
SW^) = Sn(u>)/ < A 3 > 2 

The noise terms In (su , lN {m , r) , In ( pt ^ s) are mutually independent, so that 


C / \ , c* , \ , c , a / \ . tt?5’i(a») + Q 2 5 2 (w) + a^5 , 3 (a;) 

Sgw{u) = S( 4 „,)(u>) + S( mir) {u>) + 5( prei ) + S( rp )M + — >2 

+ + Sa{u J I 5 " (u,) 


where Si,S 2 , S 3 are the Fourier transforms of the convolutions <: Ai,Aj :>,<: A 2 ,A 2 :>,<: 
8A 3 ,8A 3 :>, while Si 2 and S 2 i represent the Fourier transforms of the convolutions <: Ai,Tm* 
A 2 :> and <: T\f* A 2 ,A\ :> respectively. 

The beam d 2 , entering the port 2 of our interferometer, is generated by a degenerate parametric 
oscillator (0P0) excited by a pump beam d p , obtained by duplicating the laser beam di, entering 
Pi. In the following we will treat di(t) = e x ^,Jni{t) as a classical field (semiclassical analysis) 
whose instantaneous phase <f>(t) and intensity fluctuations 8ni(t) = n/(<) — < ni >/ will be assumed 
to be both Gaussian and mutually independent stationary processes, with autocorrelations < 
(<£(t) - <£(0)) 2 >= <tJ(|t|), < dn/(r),dn/(0) >= afCi(\r\) with of =< ( 8n t ) 2 >. 

The evolution of the the field operator d 2 has been derived by Collett and Gardiner [13] for a 
classical coherent pump. We have integrated C.-G. equation of motion of d 2 by representing the 
pump as a p = T)e 2,e+2ut> ni and applying the WKB method. 

The expectation values of the Eq.(14) have been obtained by averaging over the noise entering 
the OPO and the laser field amplitude and phase. 

In particular, as a consequence of the classsical approximation for d p we can write 


Ai = d (1) ^a (2) + d (2) td (1) = aiXjf\ A 2 - a t X^j 2 , A 3 = n t - n (2) 


(16) 


with n^ 2 l = a^ 2 \ 
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